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Climatology of Mid- latitude Ionospheric Disturbances 
from the Very Large Array Low-frequency Sky Survey 

J. F. Helmboldt/ W. M. Lane/ & W. D. Cotton^ 

The results of a climatological study of ionospheric disturbances derived from 
observations of cosmic sources from the Very Large Array (VLA) Low-frequency Sky 
Survey (VLSS) are presented. We have used the ionospheric corrections applied to the 74 
MHz interferometric data within the VLSS imaging process to obtain fluctuation spectra 
for the total electron content (TEC) gradient on spatial scales from a few to hundreds of 
kilometers and temporal scales from less than one minute to nearly an hour. The 
observations sample nearly all times of day and all seasons. They also span latitudes and 
longitudes from 28° N to 40° N and 95° W to 114°W, respectively. We have binned and 
averaged the fluctuation spectra according to time of day, season, and geomagnetic {Kp 
index) and solar (F10.7) activity. These spectra provide a detailed, multi-scale account of 
seasonal and intraday variations in ionospheric activity with wavelike structures detected 
at wavelengths between about 35 and 250 km. In some cases, trends between spectral 
power and Kp index and/or F10.7 are also apparent. In addition, the VLSS observations 
allow for measurements of the turbulent power spectrum down to periods of 40 seconds 
(scales of ~0.4 km at the height of the E'-region). While the level of turbulent activity 
does not appear to have a strong dependence on either Kp index or F10.7, it does appear 
to be more pronounced during the winter daytime, summer nighttime, and near dusk 
during the spring. 



1. Introduction 

The ionosphere is known to contain a host of tran- 
sient ionospheric disturbances. They range from flne- 
scale, turbulent like density fluctuations [e.g., Cohen 
and Rottgering, 2009] to large-scale plasma waves 
[e.g., Borries et a/., 2009] and large areas of storm- 
enhanced density [e.g., Foster et a/., 2002]. Some of 
the most common among these are so called travel- 
ing ionospheric disturbances (TIDs) which are known 
to be plasma density waves with a range of sizes 
and speeds. TIDs are often divided into two classes, 
large-scale TIDs (LSTIDs) and medium-scale TIDs 
(MSTIDs) . LSTIDs are typically associated with ge- 
omagnetic storms and have sizes ~ 2000 km, periods 
of ^1 hour, and speeds of about 700 m s~^ [Bor- 
ries et al, 2009]. MSTIDs have periods between 10 
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minutes and ^ 1 hour, speeds between 50 and a few 
hundred m s~^, and sizes ^ 100 km [e.g., Jacobson 
and Erickson, 1992; Jacobson et a/., 1995]. 

Unlike their large-scale counterparts, MSTIDs are 
not typically associated with geomagnetic storms 
and are likely generated by more localized processes. 
While the climatology of these relatively commonly 
occurring waves has been studied in the past [see, 
e.g., Hernandez- Paj ares et a/., 2006; Tsugawa et a/., 
2007], the mechanisms responsible for generating 
them are still not clear. They are most frequently 
seen at dawn and dusk, implying a relationship to 
the dynamics associated with the transition between 
the daytime and nighttime ionospheres. In the north- 
ern hemisphere, daytime MSTIDs are seen most of- 
ten in the winter, usually propagating toward the 
southeast [Hernandez- Paj ares et a/., 2006]. MSTIDs 
that appear during the night are most prevalent 
in the summer, mostly directed toward the south- 
west [Hernandez- Paj ares et a/., 2006; Tsugawa et a/., 
2007]. The nighttime MSTIDs in particular remain 
something of a mystery, especially regarding their 
orientation and direction of motion. Their formation 
may be related to the Perkins instability [Perkins, 
1973], and/or to coupling between instabilities within 
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the E and F regions [see, e.g., Cosgrove and Tsun- 
oda, 2004]. Different mechanisms have been pro- 
posed to explain their southwestward motion [Kelley 
and Makela, 2001; Cosgrove, 2007; Yokoyama et al., 
2009], however, none have been conclusively verified 
by observations. 

In addition, MSTIDs have been observed to prop- 
agate in a variety of other directions, attributed to 
changes in the F-region neutral wind and/or oro- 
graphic gravity waves [e.g., Shiokawa et a/., 2008; 
Dymond et al, 2011; Helmboldt et a/., 2012b] and to 
deposition of material from the plasmasphere [Helm- 
boldt and Interna, 2012]. Both of these types of 
MSTIDs in particular can have relatively weak am- 
plitudes, between 10"^ and 10"^ TECU (1 TECU= 
2q16 ■jrR~'^). The presence (or lack thereof) of 
MSTIDs may also be linked to increased levels of 
turbulent activity within the ionosphere [e.g., Helm- 
boldt et a/., 2012b]. 

A detailed understanding of the generation, main- 
tenance, and evolution of MSTIDs requires a clima- 
tological study which is sensitive to the full range 
of observable MSTIDs and to smaller-scale density 
fluctuations to which they may be related. Recently, 
Helmboldt et al. [2012a, b] have demonstrated that 
radio interferometers, especially those observing in 
the VHP regime, may be powerful tools for produc- 
ing such a study. With the capability to measure 
differential total electron content (TEC) with a pre- 
cision < 10~^ TECU on spacings of 1 km or less 
when observing a bright cosmic source [Cohen and 
Rottgering, 2009; Helmboldt et a/., 2012a], such in- 
terferometers are capable of observing small TEC 
fluctuations such as those related to turbulence. Us- 
ing long temporal baselines, Helmboldt et al. [2012a] 
showed that one can also detect and characterize the 
properties of MSTIDs as they pass over the interfer- 
ometer, some of which have relatively small fluctua- 
tion amplitudes (~10~^ TECU). In addition, Helm- 
boldt and Intema [2012] have demonstrated that by 
using simultaneous observations of several moder- 
ately bright VHF sources, one can characterize fluc- 
tuations on much larger scales. For instance, at 74 
MHz, the 15° field of view of the Karl G. Jansky 
Very Large Array (VLA) in New Mexco corresponds 
to an area 80 km in diameter at a height of 300 km. 

Here, we seek to apply these radio interferometer- 
based techniques to a large, 74 MHz survey con- 
ducted with the VLA to obtain a statistical, climato- 
logical picture of ionospheric disturbances spanning 



a wide range of scales in size, oscillation period, and 
amplitude. In terms of amplitude and size, this study 
constitutes one of the largest dynamic range inves- 
tigations of its kind, and it contains one of the first 
climatological studies of ionospheric turbulence. In 
Section 2, we briefly describe the characteristics of 
the survey used, the VLA Low-frequency Sky Survey 
[VLSS; Cohen et a/., 2007]. In the next section, we 
detail the analysis used and ionospheric data prod- 
ucts that were generated. In Section 4, we show the 
results of our climatological study and discuss them 
further in Section 5. 

2. Experimental Detail 

In order to create high quality astronomical im- 
ages, the VLSS standard data processing corrects 
the data for ionospheric effects. In this study, we 
re-analyze the data products from that processing to 
extract information about the dynamics of the iono- 
sphere. 

2.1. The Survey 

The VLSS is a survey of the northern sky covering 
declinations north of about —30°. It was conducted 
between 2001 and 2007 using the 74 MHz system 
of the VLA [Kassim et a/., 2007] with full dual po- 
larization (VLA program numbers AP397, AP441, 
AP452, and AP509). The VLA was used in its B- 
configuration, an inverted "Y" of 27 antennas span- 
ning roughly 11 km. Observations of sources with 
low declinations were conducted with the VLA in the 
hybrid "BnA" configuration which has the northern 
arm extended to its A-configuration size of nearly 
20 km. Hybrid configurations such as BnA are rou- 
tinely used when observing sources at low elevations 
to improve the uniformity of the antenna positions 
projected onto the sky. Approximately 20% of the 
VLSS observations were conducted in the BnA con- 
figuration. 

Other VLA configurations, which are available 
roughly four months at a time, are the A, C, and 
D configurations spanning roughly 40, 4, and 1 km, 
respectively. Observations made at 74 MHz with 
the two most compact configurations have a substan- 
tially larger effective system temperature due to the 
fact that their shorter baselines make them more sen- 
sitive to diffuse emission from the Galaxy, generally 
referred to as "sky noise." The Galactic emission 
is dominated by synchrotron emission that is much 
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brighter at lower frequencies and on relatively large 
angular scales. Consequently, the A and B config- 
urations have lower sky noise contributions due to 
the fact that their smallest baselines are too big to 
be sensitive to the large- angular- scale Galactic emis- 
sion. The B configuration was chosen for the VLSS 
because its angular resolution, 80 arcseconds, is sim- 
ilar to that of other northern sky surveys at 325 MHz 
[54 arcseconds; Rengelink et al.^ 1997] and 1400 MHz 
[45 arcseconds; Condon et al, 1998] to facilitate the 
compilation of radio- frequency spectra for detected 
sources. 

The details of the survey and associated data pro- 
cessing are discussed in detail by Cohen et al [2007] . 
Recently, improvements in imaging and calibration 
techniques as well as new algorithms for mitigating 
radio frequency interference (RFI) were combined to 
re-process all of the VLSS data. The results of this 
effort are referred to as "VLSS redux," or VLSSr, 
and the improvements over the previously used pro- 
cessing techniques are discussed by Lane et al [2012]. 

Briefiy, the VLSS mapped the sky using 523 over- 
lapping pointings and imaging the entire 15° wide, 
circular field of view of the VLA at 74 MHz at 
each pointing. Each pointing was observed for a to- 
tal of about 1.25 hours or more, broken into four 
or more 20- minute segments, or three or more 25- 
minute segments. For each pointing, these segments 
were spaced by hours or days. This amounted to a to- 
tal of 1862 separate ^20-minute observations. These 
observations cover nearly all times of day and season 
and the associated ionospheric pierce-points span a 
range in latitude and longitude, which we have illus- 
trated in Fig. 1. 

2.2. Field-based Calibration 

The effect of the ionosphere on the VLSS observa- 
tions varied significantly across the 15° field of view, 
which corresponds to a width of 80 km at an alti- 
tude of 300 km. The first-order effect of the iono- 
sphere is to shift the positions of observed sources 
by an amount proportional to the mean TEC gra- 
dient over the area spanned by the array along the 
line of sight to each source. The VLSS used a proce- 
dure referred to as "field-based" calibration [Cotton 
et a/., 2004] to account for these shifts as a function 
of time and position within the field of view. This 
calibration scheme uses "snapshot" images of rela- 
tively bright sources using 1-2 minutes of data for 
each snapshot. The positions of these sources are 



then measured within the snapshots and compared 
to a reference catalog. The original VLSS used a 
VLA-based, 1400 MHz catalog [Condon et al, 1998] 
as its reference. This produced a catalog of 74 MHz 
sources whose positions are tied to the bright sources 
contained within the 1400 MHz reference catalog. 

Substantial improvements were obtained in iono- 
spheric calibration by Lane et al. [2012] for the 
VLSSr by using the original VLSS as the reference 
catalog. In other words. Lane et al [2012] essen- 
tially produced a second iteration of the field-based 
imaging scheme where the catalog of 74 MHz sources 
generated by the first iteration was used as the refer- 
ence catalog during the second iteration. Since there 
is a variety of spectral shapes among cosmic radio 
sources, not all bright 1400 MHz sources are also 
bright at 74 MHz, and vice versa. Consequently, 
having a reference catalog of sources known to be 
bright (here, > 2.5 Jy; 1 Jy = 10"^^ W m-^ Hz-^) 
at 74 MHz made the field-based calibration process 
more robust in that a larger area could be searched 
within each snapshot for the shifted source position. 

2.3. Self-calibration 

While not used as part of the VLSS, a separate 
calibration technique referred to as "self-calibration" 
[Cornwell and Fomalont, 1999] can be used to derive 
ionospheric corrections toward a single bright cos- 
mic source if the intensity of that source dominates 
the field of view. Self-calibration uses a model for 
the brightness distribution on the sky (often a sim- 
ple point-source) to solve for antenna-based phase 
corrections. For N antennas, there are N{N — l)/2 
redundant pairs or baselines, making this an over- 
determined problem. The derived phase correc- 
tions are proportional to the difference in TEC, or 
"^TEC" , along the lines of sight of the antennas and 
an arbitrarily chosen reference antenna. Helmholdt 
et al [2012a] showed that using a bright source with 
the VLA at 74 MHz, these (5TEC measurements can 
be made to a precision of about 3 x 10~^ TECU, 
allowing for an analysis of very fine-scale, small am- 
plitude fiuctuations along the lines of sight toward a 
particular source. 

To make such fine-scale ionospheric measure- 
ments, we sought to apply self-calibration when- 
ever possible to the VLSS observations to com- 
plement the wide-field, medium-scale measurements 
provided by the field-based calibration process. To 
do this, we had to first determine how bright a source 
needed to be to accurately compute ^TEC using self- 
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calibration- derived phase corrections. We did this by 
first identifying ah sources in the VLSS catalog with 
peak intensities >6 Jy beam~^. We then used each 
source to perform phase-only self-calibration (i.e., we 
did not solve for the amplitudes of the complex an- 
tenna gains), assuming a simple point-source model. 
Since we would like to probe the smallest scales in 
space and time, we did this using the smallest time 
interval possible. For the VLSS data, this was 10- 
20 seconds. We then converted these self-calibration 
phases to STEC measurements by unwrapping the 
phases within each ~20 minute observations and de- 
trending with a linear fit. This was done separately 
for each polarization so that we could use the rms dif- 
ference between the de-trended ^TEC measurements 
from the two polarizations to assess the precision of 
the measurements. 

In Fig. 2, we have plotted ^TEC rms versus the ob- 
served peak intensity (i.e., modified by the antenna 
response) for these sources. Above a limit of about 
15 Jy beam~^, there is a noticeable anti-correlation 
between this rms and the source intensity which is 
what one would expect if the STEC measurements 
represent real ionospheric fluctuations and not noise. 
Below this limit, the data essentially form a scatter 
plot, and the self-calibration phases are likely domi- 
nated by noise caused by "confusion," i.e., the con- 
tribution from other sources in the field of view that 
were ignored in the self- calibration process. There- 
fore, for the analysis that will be described below, 
we have only used sources with observed intensities 
> 15 Jy beam~^. There are 281 sources which meet 
this criterion that were observed within 477 differ- 
ent VLSS pointings (recall that the pointings over- 
lapped). These yielded 1978 separate observations of 
the 281 bright sources. 

In addition to the 1978, 20-minute observations 
of relatively bright sources, the VLSS data also con- 
tain several five- minute observations of the extremely 
bright source Cygnus A, or "CygA." At 74 MHz, 
Cyg A is brighter than 17000 Jy, making it an ex- 
cellent source to use to determine corrections for 
instrumental effects which are relatively stable in 
time. Consequently, Cyg A was observed 289 times 
during the VLSS campaign, each time with a du- 
ration of five minutes. The extreme brightness of 
Cyg A also provides an excellent means for probing 
the smallest-scale/lowest-amplitude ionospheric fluc- 
tuations observable with the VLA [see Helmboldt et 
a/., 2012a, b]. Therefore, to increase the dynamic 



range of our study of ionospheric disturbances, we 
also applied self-calibration to all of the five-minute 
observations of Cyg A. We note that because the 
Cyg A observations were shorter, the de- trending of 
the STEC measurements made from these observa- 
tions biased them toward smaller-period oscillations. 
This makes the Cyg A data complementary to the 
data produced from the 281 bright source observa- 
tions described above, which are less sensitive to 
smaller-scale/smaller amplitude fluctuations because 
the sources used are not nearly as bright as Cyg A. 

3. Spectral Analysis 

3.1. Mult i- source Data 

The repository of source position shifts compiled 
during the processing of the VLSS provides a rich 
database of TEC gradient time series that can be an- 
alyzed and searched for evidence of wavelike activity. 
Indeed, Helmboldt and Interna [2012] showed that by 
Fourier transforming such a time series in time, and 
then in space, one may produce a three-dimensional 
(one temporal, two spatial) power spectrum "cube" 
of TEC gradient fluctuations. In practice, a spec- 
tral cube is generated for both components (north- 
south and east- west) of the TEC gradient. These two 
cubes are then added together such that the observed 
power of a wave with a TEC amplitude A and spatial 
frequency ^ will simply be (27r^74)^. These spectra 
are also normalized by the amplitude of the impulse 
response function (IRF), computed using the source 
positions projected to an ionospheric height of 300 
km. This is denoted in the units reported for these 
spectra by including a unit of IRF~-^. 

The three-dimensional spectral cubes can be used 
to map the level and direction of any detected wave- 
like structures within the data. We have applied the 
spectral analysis algorithm described in Helmboldt 
and Interna [2012] to each roughly twenty- minute ob- 
servation of each VLSS field using outputs from the 
VLSSr, yielding 1862 power-spectrum cubes. Since 
each spectrum is based on position shift data from 
multiple sources, for the remainder of this paper, we 
will referred to these spectral data as "multi-source" 
data for convenience. 

It was demonstrated by Helmboldt and Intema 
[2012] that these spectra have excellent sensitivity, 
being able to detect fluctuations with amplitudes 
< 10"^ TECU. This is in part due to the improved 
data processing/calibration techniques detailed in 
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Lane et al [2012]. They demonstrated the effect 
of these improvements on ionospheric analysis us- 
ing mean fluctuation spectra from the entire survey. 
They showed that the use of the original VLSS as 
the reference catalog as well as the new RFI miti- 
gation techniques substantially improved the sensi- 
tivity of the power spectra to both large and small 
amplitude fluctuations. Using a better reference cat- 
alog allows one to reliably search a larger area for 
each calibration source, facilitating the detection of 
large amplitude fluctuations. Reducing the noise 
in each snapshot image by implementing new RFI- 
subtracting and flagging software allowed for the de- 
tection of smaller position shifts. Specifically, the 
spatial frequency at which the mean fluctuation spec- 
trum reached the noise "floor" was increased by a fac- 
tor of two by the RFI- mitigation software, decreasing 
the typical wavelength of the weakest detectable dis- 
turbances from 70 to 35 km. 

3.2. Single-source Data 

While the power spectrum cubes described above 
provide a wealth of information about TEC fluctua- 
tions on a range of scales, their sensitivity is gener- 
ally limited at relatively high spatial and temporal 
frequencies due to the effective smoothing applied 
to the data. The position shifts of the calibrator 
sources give the mean TEC gradient over the span 
of the VLA, which in the case of the VLSS amounts 
to smoothing the data with an 11-km wide circular 
filter. In addition, snapshots were made using 1-2 
minutes of data, resulting in temporal smoothing. 

However, the measurements of (^TEC made us- 
ing individual bright sources described in Section 
2.3 used temporal sampling of 10-20 seconds with 
mean antenna spacings < 1 km. This implies that 
there is potentially more information within these 
data about fluctuations on scales too small to be de- 
tected with the multi-source spectra. To help ex- 
plore phenomena with smaller amplitudes and sizes, 
we sought to apply the methods developed by Helm- 
boldt et al. [2012a, b] and refined by Helmboldt and 
Interna [2012] to use 74 MHz VLA observations of a 
single bright source to characterize TEC fluctuations 
on significantly smaller scales. 

The single-source spectral analysis is similar to 
that of its multi-source counterpart in that a three- 
dimensional Fourier transform of the TEC gradient is 
used. However, in this case, the gradient is measured 
at each antenna for a single source rather than to- 



ward multiple sources. Because of the Y-shape of 
the VLA, the STEC measurements obtained from 
self-calibration (see Section 2.3) cannot be used to 
directly compute the full two-dimensional TEC gra- 
dient at each antenna. However, Helmboldt et al. 
[2012a] showed that most of the structure in the TEC 
surface observed by the VLA can be recovered with 
a second-order polynomial fit to STEC as a function 
of antenna position. 

Using these polynomial fits to estimate the gra- 
dient at each antenna, the same Fourier-based pro- 
cedure can be applied as was used with the multi- 
source approach. The details of how this is imple- 
mented are discussed at length by Helmboldt and In- 
terna [2012]. There are two main difference between 
the spectral cubes produced from the single source 
data and those generated using multi-source data. 
First, the spectral resolution in the spatial regime of 
the single-source-based spectra is much lower given 
the compact size of the VLA (11-km). Because of 
this, each temporal mode within each spectral cube is 
well approximated by a single spatial frequency and 
direction. Consequently, for each single-source spec- 
trum, we computed a single peak power and weighted 
(with the spectral power) mean spatial frequency in 
the north-south, ^ns, and east- west, ^ew, directions 
for each temporal mode. When combining all tem- 
poral frequencies together, this allows for a measure- 
ment of the distribution of fluctuation strength in the 
^NS^£,EW plane with much better spectral resolution 
than could normally be achieved with such a com- 
pact array. This is because, by using a single spatial 
mode for each temporal mode, we have essentially 
converted the relatively long temporal baselines to 
spatial ones. For instance, for a TID with a speed of 
100 m s~^, a 20-minute baseline amounts to a 120- 
km baseline, more than ten times larger than the 
B-configuration VLA. 

The second difference between the single- and 
multi-source analyses is that because the cosmic 
sources are essentially infinitely far away, the lines 
of sight of the antennas toward a single source are 
basically parallel. The lines of sight toward multi- 
ple sources, however, are far from parallel. Thus, 
the physical separation among pierce-points for the 
multi-source observations is proportional to the as- 
sumed height, whereas for the single-source data, the 
pierce-point separations remain virtually the same 
for any assumed altitude. This renders the multi- 
source data incapable of sensing disturbances in the 
upper ionosphere and plasmasphere that have wave- 
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lengths less than about 1000 km. Conversely, the 
single-source observations are sensitive to ^ 10- km 
and larger sized fluctuations, evan at plasmaspheric 
heights [see, e.g., Helmboldt and Interna, 2012]. 

The single-source spectral analysis was applied to 
all 20-minute observations of the 281 identifled bright 
sources and all five-minute observations of Cyg A (see 
Section 2.3). In doing so, we found that for the obser- 
vations conducted in the hybrid BnA configuration 
(see Section 2.1), the extended, 20-km long northern 
arm of the VLA posed a problem for the polyno- 
mial fit-based analysis. In short, the polynomial fits 
only provide a uniform estimate of the TEC gradient 
across the array when it is symmetric. Fitting a two- 
dimensional polynomial to an asymmetric array can 
lead to biasses within the spectral analysis performed 
using the results from such polynomial fits. For in- 
stance, we found that data from this configuration 
produced two-dimensional fluctuation spectra that 
were, on average, elongated along the north/south 
axis. Because of this inherent bias, single-source BnA 
observations (360 for the bright sources; 87 for Cyg 
A), were excluded from this analysis. 

In addition to this three-dimensional spectral anal- 
ysis, we also used the single source data to yield a 
statistical description of any isotropic fluctuations on 
the flnest scales observable within the VLSS data. 
This was done according to the prescription of Helm- 
boldt et al. [2012a] for using the self-calibration- 
determined ^TEC values to numerically compute the 
projection of the TEC gradient along each VLA 
"arm" at each antenna. While directional informa- 
tion is essentially lost within this approach, one may 
use it to recover information about isotropic fluc- 
tuations on scales < 1 km. To this end, we have 
computed a median power spectrum among all an- 
tennas from each single-source observation by sim- 
ply performing a Fourier transform of the projected 
TEC gradient time series for each antenna and deter- 
mining the median spectral power at each temporal 
frequency. Using the median rather than the mean 
helps to mitigate the influences of both strong wave- 
like disturbances and any spurious data that might 
be present. These median one-dimensional spectra 
therefore provide an excellent means to explore the 
climatological nature of turbulent fluctuations, which 
is quite complementary to the wave-based analy- 
sis provided by the multi-source spectral cubes and 
single-source, two-dimensional spectral maps. 



4. Climatology 

4.1. Mult i- source 3-D Spectra 

As described above, we used the outputs from 
the VLSSr fleld-based calibration procedure for each 
of the 1862, ^20-minute observations to construct 
fluctuation spectrum cubes, one temporal dimension 
and two spatial, of the TEC gradient. We then 
binned these by local time and then within local time 
bins, by day of the year, Kp index, and F10.7. We 
used four time of day bins, dawn (04:00-08:00), day- 
time (08:00-16:00), dusk (16:00-20:00), and night- 
time (20:00-04:00). We likewise used four day of the 
year bins, spring (34-124), summer (125-215), au- 
tumn (216-306), and winter (307-33). During the 
observations, the Kp index varied between and a 
little more than 5, and we consequently designed 
our four Kp bins to span the range 0-5. The val- 
ues of F10.7 ranged from about 70 Solar Flux Units 
(SFU; 1 SFU = 10-22 W m-^ Hz-^) to more than 
220 SFU between 2001 and 2007 (i.e., roughly so- 
lar maximum to solar minimum). We therefore con- 
structed four logarithmically spaced bins spanning 
log F10.7 = 1.85-2.35. 

In Fig. 3, we display the peak power over all tem- 
poral frequencies as a function of ^ats and £,ew of 
the average power spectrum cube for each time-of- 
day/time-of-year bin. Maps of the accompanying 
peak temporal frequency are shown in Fig. 4. Within 
each panel, the number of observations used in the 
average is printed. For all of the spectra, the full 
width at half maximum (FWHM) of the IRF is 0.009 
km"^. Within the spectral maps, there are a vari- 
ety of features that are > the size of the IRF which 
likely represent real wavelike structures rather than 
random fluctuations. For example, if one examines 
the spring nighttime panel in Fig. 3, one can see a 
prominent feature in the lower left quadrant of the 
plot. These maps are constructed such that north is 
up and east is to the right, implying that this fea- 
ture corresponds to southwestward-directed waves. 
There is a polar grid plotted as dotted grey lines to 
help one estimate direction and wavelength. In this 
case, the direction is about 45° south of west and the 
spatial frequency, ^, is about 0.015 km"^, implying a 
wavelength of 67 km (i.e., 1/0- If examines the 
same region in the corresponding panel in Fig. 4, on 
can see that the peak temporal frequency is about 5 
hr~-^. Taken with the wavelength, this gives a speed 
of roughly 330 km hr~^, or about 90 m s~^. 
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Some of the features apparent within the panels 
of Fig. 3 and 4 seem to be specific to a given season 
and/or time of day. Some coincide with the known 
climatological behavior of previously explored fluc- 
tuations. A summary of prominent features found 
within these data and the single-source spectral anal- 
ysis presented below is given in Table 1, including 
comparisons with previous work. 

As alluded to in Section 1, southwestward-directed 
MSTIDs have been previously observed to be promi- 
nent during summer nighttime in North America. 
Likewise, winter daytime MSTIDs are known to be 
common in the same region, usually propagating to- 
ward the southeast [e.g., Hernandez- Paj ares et a/., 
2006; Tsugawa et a/., 2007]. We see direct evi- 
dence of these same phenomena within the multi- 
source spectra. The summer nighttime spectrum 
shows a southwest-directed feature with a wave- 
length of about 100 km and a period of about 17 
minutes, implying a speed of 100 m s~-^. This 
is roughly consistent with the known properties of 
MSTIDs. Similarly, the winter daytime spectrum 
shows a southeastward-directed feature with a peak 
at a wavelength of 190 km, a period of 17 minutes, 
and an implied speed of 190 m s~^, which is again, 
consistent with MSTIDs. 

There are also many apparent features within the 
multi-source spectra that are not commonly known. 
For instance, several show evidence of wave activity 
directed toward the northeast, namely spring dawn, 
summer daytime and dusk, and autumn dusk and 
nighttime. There are southwestward-directed waves 
during spring nighttime, but they are smaller (^50 
km) and slower (^60 m s~^) than the MSTIDs seen 
during summer nighttime. In addition, the spring 
nighttime spectrum shows evidence of MSTID-like, 
westward- directed waves not seen during other time 
periods. While medium-scale wave activity seems to 
be common during the daytime, the predominant di- 
rection appears to have a seasonal dependence with 
southwestward-directed waves being more prevalent 
during summer daytime and southeastward-directed 
waves dominating during autumn and winter. The 
autumn dusk and nighttime spectra are particu- 
larly interesting as they show a plethora of features 
with no apparent dominant direction or size/speed. 
Most of the spectra show evidence of relatively small 
(1^1 > 0.01 km~^) features in the periphery of the 
areas they occupy within the £,ns,£,ew plane. 



In Fig. 5, we have plotted spectra similar to those 
displayed in Fig. 3 but for bins of time-of-day and 
Kp index. Once can see from these that most of 
the features seen in Fig. 3 do not appear to have a 
clear dependence on geomagnetic activity. For in- 
stance, southeastward-directed waves are seen dur- 
ing the daytime at all values of Kp. However, it does 
appear that the northeastward-directed waves that 
appear at dawn are more prominent at low to mod- 
erate values of Kp as are the southwestward-directed 
nighttime waves (both the ~50 and ^100 km sized 
waves) . 

Similar spectra are plotted again in Fig. 6, this 
time for bins of time-of-day and log F10.7. Like 
the results for the Kp index, most of the features 
seen in Fig. 3 do not seem to have noticeable depen- 
dences on solar activity. Again, the most notable ex- 
ceptions are the dawn northeastward-directed waves 
and nighttime southwestward-directed waves, both of 
which are seen most prominently at the lowest levels 
of solar activity. The larger southwestward-directed 
waves seem to actually be most prominent at the 
highest solar activity levels. However, we note that 
all of the summer observations were taken in 2002 
near solar maximum. Indeed, the spectrum in the 
high log F10.7, nighttime bin looks very similar to 
the summer nighttime spectrum in Fig. 3. The fact 
that these waves have been seen by others to have 
a strong seasonal dependence suggests that it is the 
time of year rather than solar activity levels that has 
had the strongest influence on the appearance of this 
feature in this case. 

For the lowest Kp and log F10.7 bins, the night- 
time spectra show a significant eastward-directed fea- 
ture that is not prominent in the nighttime spectra 
in Fig. 3. The feature is actually closer to magnetic 
east, which is about 12° south of due east at the 
VLA. The wavelength and speed of this feature are 
about 150 km and 150 m s~^, respectively. Because 
these waves are predominantly magnetic eastward- 
directed with relatively large speeds, they may in 
fact be field-aligned irregularities within the plas- 
masphere similar to those discovered with the VLA 
by Jacobson and Erickson [1992]. For instance, co- 
rotation within the plasmasphere will yield an ob- 
served speed for such irregularities of 150 m s~^ for 
heights of about 2000 km. For these observations, 
this corresponds to a Mcllwain L-parameter of about 
2.1, which is typical for the irregularities previously 
discovered with the VLA [Hoogeveen and Jacobson^ 
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1997]. These disturbances will be discussed further 
in Section 5. 

4.2. Single-source 2-D Spectra 

As discussed in Section 3.2, the spectral analy- 
sis of each single-source observation yielded a mea- 
surement of the distribution of spectral power in the 
£,NS,£,EW plane that reaches higher spatial frequen- 
cies than the multi-source spectra. We have binned 
these spectral maps according to local time, season, 
Kp index and F10.7 in the same way as the multi- 
source spectra and have displayed the results in Fig. 
7-9. In each panel of each of these figures, we have 
plotted the mean power over all temporal frequencies 
from the corresponding mean multi-source spectral 
cube as blue contours. 

The spectral maps shown in Fig. 7 show that at 
all times of day and year, there are significant fiuc- 
tuations at scales too small to detect with the multi- 
source data. However, when averaged together, they 
seem to be largely isotropic with a few notable ex- 
ceptions. First, the map for winter daytime seems to 
be especially asymmetric, showing an over-density in 
the northwest quadrant. This northwest over-density 
is dominant on scales of about 100 km. Fig. 8 and 
9 show that the northwest over-density seen during 
winter daytime is more prominent at moderate solar 
activity levels, but does not seem to have a depen- 
dence on Kp. 

The autumn dawn map shows a population 
of nearly northward-propagating disturbances with 
wavelengths of about 50 km not seen in the corre- 
sponding multi-source spectrum. There is some in- 
dication from Fig. 9 that these waves occur more fre- 
quently during times of relatively high solar activity. 

Perhaps the most striking feature is a significant 
group of small-scale waves directed toward the south- 
east during spring dawn with wavelengths of about 
40 km which represent a distinct population of waves 
not detectable with the multi-source data. Fig. 8 and 
9 show that these waves also appear prominently at 
times of low geomagnetic and solar activity. 

We also discussed in Sections 2.3 and 3.2 sepa- 
rate five-minute calibration observations that were 
made of the extremely bright source, Cyg A, which 
we also used to produce spectral maps. There were 
significantly fewer of these observations (202 versus 
1618 in the B configuration), and their time-of-day 
and seasonal coverage is not as good. However, as 
also discussed in Section 3.2, the spectral maps pro- 



duced from these observations are more sensitive to 
smaller-scale fluctuations, making them complemen- 
tary to the single-source maps shown in Fig. 7-9. 
Therefore, we have displayed in Fig. 10-12 Cyg A- 
based spectral maps binned by local time and season 
(Fig. 10), Kp index (Fig. 11), and F10.7 (Fig. 12). 

The Cyg A-based spectra show no signs of any 
signiflcant groups of waves. However, as expected, 
they typically extend to larger spatial frequencies 
than their counterparts shown in Fig. 7-9. While 
the maps often show roughly isotropic distributions 
of fluctuations, reminiscent of turbulence, many show 
some evidence of "preferred" directions. For in- 
stance, the spring nighttime map is elongated along 
the northwest-southeast axis with more power appar- 
ent in the southeast direction. The same appears to 
be true during dawn for higher Kp indices, but with 
a much larger fraction of the total spectral power 
in the southeastern quadrant. The summer dawn 
spectrum shows a signiflcant northwestern extension 
toward higher spatial frequencies. 

The summer dawn map has an extension toward 
the northwest that reaches spatial frequencies as high 
as about 0.04 km~^ (or, scales of 25 km). Similar 
features can be seen in the dawn maps in Fig. 12 for 
moderate to high levels of solar activity. There is also 
a prominent group of northwestward-directed waves 
extending to large (0.05 km~^) spatial frequencies in 
the spring dusk map, but this is based on a single 
observation. 

4.3. Single-source 1-D Spectra and Turbulence 

As described in Section 3.2, we have produced me- 
dian TEC gradient fluctuation spectra using the arm- 
based method of computing the projected TEC gra- 
dient from all single-source data. Since the projected 
gradients do not rely on polynomial flts, we have used 
all 1978 single-source observations for this analysis, 
including those conducted in the hybrid BnA conflg- 
uration. These one- dimensional spectra probe the 
shape of the spectrum of isotropic of fluctuations 
down to the smallest periods possible, 40 seconds. 
If we assume that the population of these fluctua- 
tions does not change signiflcantly during each ob- 
servation, the temporal frequencies can be directly 
related to spatial ones using the apparent speed of 
the observed cosmic source at ionospheric heights. 
Since turbulence is likely to be much more promi- 
nent at lower altitudes where ion-neutral coupling is 
more important, we have computed the median ap- 
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parent source speed at the height of the £^-region, 
100 km, which is roughly 10 m s~^. This imphes 
that these spectra probe scales as sraall as 0.4 kni. 

Continuing with the assumption of "frozen" turbu- 
lence, we can characterize the level of turbulent activ- 
ity by fitting a simple power-law model to each spec- 
trum. Specifically, following the work of Kolmogorov 
[1941a, b] and Tatarski [1961], the turbulent spec- 
trum of phase fluctuations, which in this case is the 
same as TEC fluctuations, should be proportional 
to Therefore, the TEC gradient spectrum 

should be cx and the temporal spectra should 

be well fit by a spectrum oc for frozen tur- 

bulence. We have displayed mean one-dimensional 
single-source spectra in Fig. 13, binned by season 
and local time, to illustrate that this is generally 
the case. For this analysis, we have excluded the 
five-minute Cyg-A observations because the shorter 
duration yields lower temporal spectral resolution 
which significantly alters the shapes of the spectra at 
lower frequencies. We also found that the 20-minute, 
single-source observations were more than adequate 
for characterizing the spectrum of turbulent fluctua- 
tions (see Fig. 13) and that the Cyg A-based spectra 
did not significantly improve or add to this analysis. 

For each single-source spectrum, we have fit a sim- 
ple model of P^i/"^/^ + where Pt is a parameter 
that characterizes the level of turbulent activity and 
TV is a constant that models the effect of noise within 
the spectrum. The fits were constrained to z/> 20 
hr~^ because below this limit, the influence of the 
window function used (i.e., a simple 20-minute wide 
boxcar) as well as relatively strong, medium-to-large 
scale fluctuations on the shape of the spectrum can 
be significant. The fits are plotted as red curves in 
Fig. 13. One can see that these turbulence model 
fits provide a good approximation of the data in all 
cases. 

To obtain a climatological picture of turbulent ac- 
tivity, we also fit the turbulence model described 
above to the arm-based spectrum from each single- 
source observation. Among these spectra, the 
parameter Pt varied between roughly 1 and 20 
(mTECU km-i IRF-i)^ and the median value for 
the parameter N was 10"^ (mTECU km"^ IRF-i)^. 
The mean value of Pt was then computed within 
bins of local time and day of the year, Kp index, and 
log F10.7. These two-dimensional distributions are 
displayed in the panels of Fig. 14. During the day, 
the turbulent activity appears to be highest during 
winter months, which is when daytime MSTID ac- 



tivity also peaks. During dusk, it seems to be higher 
in the late winter/early spring while post-midnight 
and before dawn, there seems to be some increased 
activity during the summer. There also seems to 
be a "spike" of turbulent activity during the day at 
moderately high values of Kp (4+ to 5—). Daytime 
turbulent activity seems to be relatively high during 
moderate levels of solar activity (F10.7 ^^100 SFU). 
Post-midnight and before dawn, the turbulent activ- 
ity also has peaks near FIO.7^120 and 170 SFU. 

5. Discussion 

We have used a large, 74 MHz survey of the 
northern sky, the VLSS, to explore the rich envi- 
ronment of ionospheric disturbances over the south- 
western United States and northern Mexico (see Fig. 
1). Some of the results were expected from previ- 
ous analysis, namely the prominent MSTIDs dur- 
ing summer nighttime and winter daytime, which 
propagate toward the southwest and southeast, re- 
spectively [see, e.g., Herndndez-Pajares et a/., 2006; 
Tsugawa et a/., 2007]. However, there are several fea- 
tures unique to this study, and we discuss a subset 
of them below. 

5.1. Muti-directional Waves 

The results shown in Fig. 3 show evidence of 
waves moving in several different directions on dif- 
ferent scales. Most notable are the results for au- 
tumn dusk and nighttime, each showing this variety 
within a single mean spectrum. While several fac- 
tors may influence the appearance of these waves, 
it seems most likely that they are related to oro- 
graphic gravity waves [see, e.g., Vadas and Crowley, 
2010]. This is because beneath the region of the iono- 
sphere probed by the VLSS observations, the terrain 
is quite mountainous. The is illustrated within the 
bottom panel of Fig. 1 where one can see from the 
displayed relief map that the latitudes and longitudes 
of the ionospheric pierce-points for the VLSS obser- 
vations cover the most mountainous regions of north- 
ern Mexico, New Mexico, and Arizona. In addition, 
the Rocky Mountains lie largely to the north and 
west of these observations, implying that features like 
the southeastward-directed, 40-km sized waves seen 
during spring dawn may be related to waves gener- 
ated by airflow over the Rockies. 
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5.2. Northeastward Waves 

A prominent and repeating feature within the 
spectra shown in Fig. 3 is a class of northeastward- 
directed waves. These are similar in size, and often 
in strength, to MSTIDs, but are traveling in a rel- 
atively abnormal direction for these types of distur- 
bances. Similar waves were observed with the VLA 
during nighttime observations in August, 2003 by 
Helmholdt et al [2012b]. The rehef map in Fig. 1 
does indicate that there are mountains to the south- 
west of many of the ionospheric pierce-points, im- 
plying that gravity waves may be involved. How- 
ever, these northeastward-directed waves are often 
relatively strong compared to other spectral features 
seen in Fig. 3 (see spring dawn and summer daytime 
for examples). The mountains to the southwest are 
not particularly large, especially when compared to 
the Rockies to the northwest, implying that there 
may be more at work here then vertically propagat- 
ing gravity waves. The northeastward-directed waves 
detected by Helmholdt et al. [2012b] were seen to co- 
incident with drops in F-region height measured with 
relatively nearby ionosondes data. This is similar to 
what was observed over Japan by Shiokawa et al. 
[2008]. It is therefore plausible that F-region com- 
pressions are related to the appearance of these fea- 
tures within the VLSS data. However, a thorough in- 
vestigation using contemporaneous ionosondes data 
is required to resolve this and is currently underway. 

5.3. Magnetic Eastward Waves 

As noted in Section 3.1, there are relatively 
unique, eastward-directed features in the nighttime 
spectra for the lowest Kp and log F10.7 bins shown 
in Fig. 5 and 6. These could be similar to the mag- 
netic eastward-directed waves discovered with the 
VLA and established by Hoogeveen and Jacobson 
[1997] to be plasmaspheric waves. They could also 
be associated with disturbances within plasma flow 
from the plasmasphere to the ionosphere discovered 
with the VLA by Helmholdt and Intema [2012]. The 
latter seems particularly likely because (1) they are 
only seen during the night when such flows are most 
likely to occur and (2) the spectra also show evi- 
dence of westward-directed, presumably ionospheric 
waves which were also found by Helmholdt and In- 
tema [2012] to occur simultaneously with the plas- 
maspheric disturbances. The fact that they are most 
prominent during periods of low geomagnetic and so- 



lar activity may imply that a relatively undisturbed 
plasmasphere is required to produce these features. 
A more detailed study of this phenomenon is cur- 
rently being conducted using GPS-based TEC maps, 
Kp and AE indices, and VLSS data. The prelimi- 
nary results suggest that these disturbed flows may 
be triggered by localized depletions within the night- 
time ionosphere triggered by forcing from the lower 
atmosphere, likely related to tides. The full results 
will be presented in a subsequent paper. 

5.4. Turbulence 

The results shown in Fig. 13 confirm what was 
demonstrated by Cohen and Rottgering [2009] and 
[Helmholdt et a/., 2012b], that the median behavior 
of the spectrum of ionospheric fiuctuations is turbu- 
lent. Fig. 14 shows that while the level of turbulence 
does vary over time, the difference between the max- 
imum and minimum levels is only about a factor of 
five. Contrast this with the results for wave activity 
displayed in Fig. 3-12 where the range in spectral 
power spans at least two orders of magnitude. The 
level of turbulent activity does appear to be related 
in some way to the occurrence of certain wave phe- 
nomena. In particular, Pt is relatively large when 
MSTIDs are prominent in winter daytime and sum- 
mer nighttime. 

6. Conclusions 

By using a large database of 74 MHz observations 
of cosmic sources, we have been able to probe the cli- 
matological behavior of a variety of ionospheric dis- 
turbances over the southwestern United States and 
parts of Mexico. We have shown that small-scale, 
turbulent fluctuations are present nearly all of the 
time with only a weak dependence on time of day, 
time of year, and/or the level of geomagnetic and so- 
lar activity. We found a variety of small (roughly 40 
km wavelength) and medium (100 km and larger) 
scale phenomena present at different times. The 
sources of these disturbances may range from gravity 
waves, plasmaspheric interactions, coupling between 
different ionospheric layers, and plasma instabilities 
(e.g., the Perkins instability). Several of these wave- 
like irregularities have relatively small amplitudes 
too weak to detect with other methods such as GPS- 
based TEC measurements. In some cases, the wave 
activity was coincident with increases in turbulent 
activity, namely during winter daytime and summer 
nighttime. 
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We note that in general, the impact of the results 
presented is somewhat hampered by the limited sea- 
sonal and time-of-day coverage of the VLSS shown 
in the upper panel of Fig. 1. For instance, all of the 
summer observations were conducted in the month of 
June. The fact that the summer nighttime MSTIDs 
appear weaker than their winter daytime counter- 
parts may be heavily influenced by the lack of data 
for July and August. In addition, the only observa- 
tions conducted during the dawn hours in the sum- 
mer were calibration observations of Cyg A. We are 
now designing a survey with the new 330 MHz system 
for the VLA, which is currently being commissioned. 
This new survey will be scheduled to optimize both 
time-of-day and seasonal coverage. It will benefit 
from two new arrays of GPS receivers operating con- 
tinuously within New Mexico which have recently in- 
creased the number of such stations from ~6 during 
the time of the VLSS to nearly 40. There will also be 
a new digital ionosonde operating in nearby Kirtland 
Air Force Base in Albuquerque which will potentially 
provide information about the height where detected 
disturbances occur. Thus, we will have a contempo- 
raneous GPS-based account of the largest fluctua- 
tions, providing a complete and unique inventory of 
the seasonal dependence of ionospheric disturbances 
that will expand upon the work presented here. 
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Table 1. Summary of Observed Phenomena 



Phenomenon Season(s) Time(s) of Day Notes 



SW-directed MSTIDs 


summer 


daytime-nighttime 


previously detected with GPS and airglow imagers; e.g 
Hernandez- Paj ares et al. [2006]; Tsugawa et al. [2007] 


SE-directed MSTIDs 


autumn/ winter 


daytime 


same as above 


NE-directed MSTIDs 


spring/summer 


dawn / daytime-dusk 


isolated cases detected previously by 
Shiokawa et al [2008] and Helmholdt et al [2012b] 


multi-directional, 
medium-scale waves 


autumn 


dusk-nighttime 


likely associated with gravity waves 


SW- and W-directed 
medium- scale waves 


spring 


nighttime 


seen with magnetic east-directed waves at low Kp] 
isolated case detected by Helmholdt and Interna [2012 


SE-directed 40-km- 
scale waves 


spring 


dawn 


found when Kp is low; 
direction close to magnetic east 


turbulence 


spring/ summer/ 
winter 


dusk / nighttime / 
daytime 


expands on diurnal variation 
seen by Cohen and Rottgering [2009] 
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Figure 1. Upper: The local time and day of the 
year of each of the VLSS observations; points are 
color-coded by year. Lower: The latitude and longitude 
of the ionospheric pierce-point for the center of the 
field of view for each VLSS observation assuming a 
height of 300 km. The same color-coding by year 
is used as in the upper panel. A relief map from 
http:/ /www. ngdc.noaa.gov/mgg/global/global. html is 
also displayed for reference. 




Figure 2. The rms difference in 6TEC between the two 
polarizations among all baselines and time steps for each 
observation of each VLSS source brighter than 6 Jy versus 
its observed peak intensity (i.e., modified by the antenna 
response). The vertical grey line shows the chosen limit 
of an intensity of > 15 Jy beam~^ for single-source data. 
See Section 2.2 for more discussion. 
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Figure 3. From mean fluctuation spectrum cubes, de- 
rived from multi-source data and binned by local time 
and season (see Section 3.1), the peak power over all 
temporal frequencies as a function of north-south and 
east- west spatial frequencies. Locations where the mean 
power over all temporal frequencies is lower than 5 x MAD 
are set to zero, where MAD is the median absolute de- 
viation among the mean power values of all spatial fre- 
quencies. The number of observations averaged together 
to make each mean spectrum is printed in each panel. 
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Figure 4. From mean fluctuation spectrum cubes, de- 
rived from multi-source data and binned by local time 
and season (see Section 3.1), the peak temporal frequency 
as a function of north-south and east-west spatial fre- 
quencies. Locations where the mean power over all tem- 
poral frequencies is lower than 5 x MAD are set to zero, 
where MAD is the median absolute deviation among the 
mean power values of all spatial frequencies. The num- 
ber of observations averaged together to make each mean 
spectrum is printed in each panel. 
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Figure 5. The same as Fig. 3, but for spectra binned by local time and Kp index. 
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Figure 6. The same as Fig. 3, but for spectra binned by local time and log F10.7. 
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Figure 7. Mean maps of TEC gradient fluctuation 
power, derived from single-source data and binned by lo- 
cal time and season (see Section 3.2). The blue contours 
represent the mean power over all temporal frequencies 
from the corresponding spectral cubes shown in Fig. 3 
and 4. The number of observations averaged together to 
make each mean map is printed in each panel. 
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Figure 8. The same as Fig. 7, but for spectra binned by local time and Kp index. 
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Figure 9. The same as Fig. 7, but for spectra binned by local time and log F10.7. 
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Figure 10. Mean maps of TEC gradient fluctua- 
tion power, derived from five-minute observations of the 
bright cahbration source, Cygnus A. These maps were 
made by binning the observations by local time and sea- 
son (see Section 3.2). The magenta contours represent 
the mean power over all temporal frequencies from the 
corresponding spectral cubes shown in Fig. 3 and 4. The 
number of observations averaged together to make each 
mean map is printed in each panel. 
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Figure 11. The same as Fig. 10, but for spectra binned by local time and Kp index. 
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Figure 12. The same as Fig. 10, but for spectra binned by local time and log F10.7. 
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Figure 13. Within bins of local time and season, the 
mean one-dimensional fluctuation spectra derived from 
single-source data using the arm-based method described 
in Section 2.2. The turbulence model described in Section 
3.3 was fit to each of these spectra and is plotted as a red 
curve. 
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Figure 14. The strength of the turbulent fluctuation 
power, Pt (see Section 3.3), derived from single-source 
data as a function of local time and day of the year (up- 
per), Kp index (middle), and F10.7 (lower). 



